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The dc Josephson effect is investigated in a single-walled metallic carbon nanotube connected to 
two superconducting leads. In particular, by using the Luttinger liquid theory, we analyze the effects 
of the electron-electron interaction on the supercurrent. We find that in the long junction limit the 
strong electronic correlations of the nanotube, together with its peculiar band structure, induce 
oscillations in the critical current as a function of the junction length and/or the nanotube electron 
filling. These oscillations represent a signature of the Luttinger liquid physics of the nanotube, for 
they are absent if the interaction is vanishing. We show that this effect can be exploited to reverse 
the sign of the supercurrent, realizing a tunable vr -junction. 



I. INTRODUCTION 

Since the discovery by lijima in 1991—, carbon nan- 
otubes have attracted much interest in the community 
of Mesoscopic Physics. Due to their peculiar electronic 
and mechanical properties, they are regarded as optimal 
candidates for nanotechnological implementations, and 
have been successfully applied to the realization of 
quantum transistors^i^, electron waveguides*, inter- 
ferometric devices^i^ as well as nanoelectromechanical 
systems^. Recent experiments have spurred the interest 
in superconducting properties of these materials: it has 
been observed indeed that proximity-induced supercon- 
ductivity can arise in nanotube bundles contacted to 
superconductors (S); in ropes, intrinsic superconductiv- 
ity has also been measured^i^ and explained in terms 
of combination of electron coupling to the breathing 
phonon modes and intertube Cooper-pair tunneling^. 
Individual multiwall nanotubes have recently been 
utilized in the fabrication of superconductor-nanotube- 
superconductor hybrid structures, allowing to reveal 
the dynamics of multiple Andreev reflectionsiS and 
to realize a controllable supercurrent transistor—. By 
contrast, the investigation of superconducting properties 
of single-walled nanotubes in hybrid structures has been 
only partly explored so far. 

Metallic Single-walled carbon nanotubes (SWNT) are 
known to behave as one-dimensional (ID) conductors 
with four conduction channels exhibiting ballistic trans- 




Figure 1: Schematic set-up of the S-TSWNT-I-S junction un- 
der investigation. 



port up to several /zmiiii^. Differently from other ID 
metals, SWNT preserve their conduction properties even 
at very low temperature, since the cylindrical lattice 
geometry prevents the arising of Peierls distorsion. 
They thus offer promising features for interconnecting 
components of nanodevices. Due to their ID char- 
acter, electronic correlations have dramatic effects on 
the behavior of SWNT: experimental evidences of a 
power-law behavior for the conductance as a function 
of temperatureii indicate that SWNT exhibit a Lut- 
tinger liquid (LL) behavior, and that their elementary 
excitations are not fermionic quasi-particles like in 
normal 3D metals^^ii^. It is thus expected that, when a 
SWNT is contacted to S leads at equilibrium, electronic 
correlations might significantly modify the behavior of 
the supercurrent with respect to junctions realized with 
a normal metal. This issue has been addressed in the 
literatur o-^^d''dSd9, 20^21 ^^^^ t^^^ been shown that the 

effect of interaction is particularly enhanced when the 
coupling between the LL and the S leads is realized 
through tunnel junctions. However, most of these inves- 
tigations focused on the case of a two-channel (i.e. one 
spinful mode) LL, and cannot be straightforwardly ap- 
plied to the case of a four-channel SWNT. In this paper 
we discuss this problem investigating the dc Josephson 
effect in a S-TSWNT-I-S junction, and show that new 
features arise due to the peculiar band structure of 
SWNT. The paper is organized as follows: in Sec. II 
we briefly review the model used to describe SWNT, 
accounting for electron-electron interaction within the 
Luttinger Liquid theory. In Sec. Ill we present our 
results about the Josephson current. We find that the 
interaction yields a twofold effect on the critical current 
jc- on the one hand it modifies the scaling law of jc as 
a function of the junction length d; on the other hand, 
it introduces oscillations of function of either 

the electron filling or the junction length d. The latter 
oscillations are absent for a non-interacting system, 
and therefore represent a signature of Luttinger liquid 
behavior on the supercurrent. Finally, in Sec. IV we 
discuss the results and propose possible implementations 
to observe this effect. 
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Figure 2: The electron band dispersion relation of a SWNT 
originates from the two-sublattices honeycomb carbon struc- 
ture, and is characterized by four Fermi points. The latter 
can be identified through two Fermi momenta: kp denotes 
the band crossing points, whereas qp accounts for the devia- 
tion from kp, i.e. the electron filling of the SWNT. 



II. MODELING THE SYSTEM 

The set-up of the system is depicted in Fig.[l] a metal- 
lic SWNT is coupled through tunnel contacts to two su- 
perconducting leads to realize a S-I-SWNT-I-S junction. 
For simplicity, we limit here our treatment to the case 
of armchair nanotubes; we also assume that the S leads 
have the same gap |A|; the two superconducting order 
parameters thus read Ai^2 = |A|e*-^i'^, where Xi is the 
superconducting phase of the «-th lead. We are inter- 
ested in the dependence of the critical current on the 
junction length d; we thus analyze the regime 



Ac , C < < i 



(1) 



where Ac represents the width of the contacts, ^ the 
coherence length of the S electrodes, d the electrode dis- 
tance, and L the length of the nanotubc. The regime ^ 
is quite realistic in view of customary fabrication of /im 
long ballistic nanotubes"'^^, and the recent realization of 
superconducting tips for Scanning Tunneling Microscope 
(STM)^i^i2i or of 10-20 nm short superconducting 
finger leads. In order to simplify the mathematical treat- 
ment without altering the essential physical features 
of the regime ([IJ, we shall henceforth assume that the 
tunnel contacts are point-like, the coherence length ^ 
is vanishing, and the length of the nanotube is infinite, 
L — !■ oo. 

In a metallic nanotube the lowest band consists of four 
electron branches located around two Fermi points akp^ 
with a = ±1; the energy separation to the second band 
is of the order of eV, so that the latter can be in prac- 
tice neglected up to a broad range of thermal excita- 
tions. Within this energy scale, the energy dispersion 
of the lower band is linear under quite good approxima- 
tion, as shown in Fig. [21 SWNT can thus be regarded 



as four-channel ID metals. As discussed in the intro- 
duction, their ID character implies that a careful treat- 
ment of the electron-electron interaction is needed. It is 
indeed well known that transport properties of SWNT 
cannot be explained in terms of the customary Fermi liq- 
uid theory, since their elementary excitations are bosonic 
plasmon modes, rather than fermionic quasi-particles. A 
model for SWNT based on the Luttinger liquid theory 
has been formulated a decade agc^^'^^, and applied in 
a number of problem a^''i^^'^^i'^°i^^i'^^'^'^ . Here we briefly 
remind the main aspects that are relevant to our discus- 
sion: A SWNT can be ideally obtained by wrapping into 
a cylinder shape a graphene sheet, whose honeycomb car- 
bon lattice consists of two sublattices p — ±. A nearest- 
neighbor tight-binding calculation of the 7r-electrons in 
the graphite, together with appropriate wrapping bound- 
ary conditions, leads to express the electron field in the 
nanotube as 



a=±,r=R/L P=± 



where a =t,i denotes the spin component and x the 
longitudinal coordinate in the nanotube. In Eq. (2]), Upr 
are the entries of the matrix 



U 



1 1 

i —i 



(3) 



describing the unitary transformation from the sublattice 
electron fields into the right (left) moving fields descrip- 
tion. The exponential terms in Eq. ^ represent the 
fast oscillating contribution to the electron wave func- 
tion, where the wave vector qp is related to the electron 
filling exceeding the Fermi points ikp, as illustrated in 
Fig. [51 Finally the field iparaix) varies slowly over the 
scale of Fermi wavelength. 

In order to account for the electron-electron interac- 
tion, it is useful to represent the electron fields ipara 
through the bosonization identity 



i^arcrix) = ^^I!- exp{z(/?aro- (a:) } 



(4) 



where iparaix) is the plasmonic field describing the 
long wavelength fluctuations. The operators rjara are 
Klein factors obeying a Clifford algebra and ensuring 
the correct anticommutation between different fermionic 
species. Finally a is a cut-off length regularizing the the- 
ory, and is of the order of the lattice spacing. The effec- 
tive hamiltonian for the SWNT reads 



/+OC 
dx 

30 



]0 



(5) 



with vp 8 ■ 10^ ms~^. Here Ojs{x) are four indepen- 
dent bosonic fields, with j labehng charge(c) and spin(s) 
degrees of freedom, and 6 = ± denoting two indepen- 
dent linear combinations of the a = ± branches. The 
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fields 4>js{x) are related to 9js{x) by the duality relation 
[0js{x),dy(pjs{y)] — iS{x — y), and Kjs are interaction 
parameters, with Kjs < 1 {Kjs > 1) for repulsive (at- 
tractive) interaction and Kjs = 1 for vanishing interac- 
tion. It can be shown'^'^ that the mode (j = c, i5 = 
is strongly interacting {K = Kc+ — 0.3), while the three 
other modes are neutral KjS:^c+ = 1- The fields faro- 
appearing in the bosonization identity ^ are linear com- 
binations of the decoupled modes 



^ara- =— {0c+ + r9c+ + aq 



(6) 



The Hamiltonian modehng the S-I-SWNT-I-S junction 
thus reads 



H ~ TisWNT + ^SCi + ^SC2 + 



(7) 



where TYswnt is given by (O, Hsci,2 ^'^'^ the usual BCS 
hamiltonians for the electrodes, and TLt describes the 
nanotube-electrodes electron tunneling. Denoting by Xi 
the nanotube coordinate of the injection point to the i-ih 
electrode, one can write 

Wt = E E r,(sl(a;,)*.(^») + *U^.)S.(a;.)) (8) 

i=l,2cr = T,i 

where sj is the electron field operator in the i-th lead, 
and Ti.2 are tunneling amplitudes. The Josephson cur- 
rent is computed perturbatively in the tunneling ampli- 
tudes. 



III. DC JOSEPHSON CURRENT 

Denoting by F the free energy of the S-I-SWNT-I-S 
junction, and by x = X2 ~ Xi the phase difference be- 
tween the two superconductors, the dc Josephson current 
is obtained as 



2edF 
h dx 



(9) 



Evaluating Ij to the fourth order in the tunneling am- 
plitudes Ti, one obtains (up to independent terms) 



F = -2 



3? 



n / 



(10) 



X Ti{ti - T2)C?(Ti,r2,T3,T4;(i)J^|(T3 - Ti)] 

In Eq. (|10p [3 denotes the inverse temperature and 



(r-r') = (T{4^(x,r)4,(x,r')}) 



7riV(0) 



lAI 



(11) 



the anomalous BCS T-ordered correlator in the i-th S 
lead, with a density of states of the normal state N{0) 
at the Fermi energy, and Matsubara frequencies uin — 
{2n + l)n/l3. Finally 



S(n,T2,T3,T4;d) 



=(T { vl^^ (0, n)*^ (0, r2)vl/| (d, ra)*! (d, T4)}) 



(12) 



is the two-electron T-ordered correlator in the SWNT. 
Under the condition ((1]), one has A » hvp/d, implying 
that the lead-nanotube tunneling time is much shorter 
than the traversal time along the junction. Eq. (|11[) 
is then well approximated by a 5(t — r'), and tunnel- 
ing effectively involves electron pairs. In this regime the 
Josephson current can be written as 



(13) 



Here the first term Iq = {2e/ h){hvF/d)T sin x accounts 
for the dependence on the superconducting phase differ- 
ence X: and corresponds to the Josephson current of a 
long ballistic junction with bare transmission coefficient 
r = (4/27r)|TiT2p7r2(jv(0)/nup)2 at zero temperature. 
The second term represents the (dimensionless) critical 
current and encodes the details of the junction: it de- 
pends on the length, on the temperature, and on the 
interaction effects, as we shall see below. Explicitly, it 
reads 



1 /dV r^/® 



e|e|)x 



(14) 



r— ±,aa' 

with 9 = kBTd/hvp- Two types of processes, denoted by 
P and A, contribute to jc- the former (latter) describes 
tunneling of Cooper pairs formed by electrons with par- 
allel (antiparallel) momenta. The related pair operators 



(15) 



Or^aa'ix,T) (a;, r)V'„,^^ (x, t) 

yield the two correlators 
C,3(fcW,6-a)- 

K2™)^(r{o,^/t(o,^a)o-/-(.,^6)}) 

(16) 

appearing in Eq. (|14p . Importantly, these two terms cor- 
respond to different dependences on the momenta defin- 
ing the four Fermi points of the SWNT: while A pro- 
cesses only involve kp, P processes are also character- 
ized by the electron filling momentum qp, as can be seen 

P I A 

from the phase factors multiplying C^.^^,, . Since typically 
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qp ^ kp, two periods are expected to arise in the depen- 
dence of the Josephson current on the junction length d. 
However, this is not necessarily the case. In the first 
instance, indeed, the dependence on kp amounts to a 
prefactor 1 + cos{2kpd), and is extremely difficult to be 
observed in a realistic system where the approximation 
of point-like contacts is not valid, for the period of these 
oscillations is usually smaller than the typical contact 
width Ac- Even in the regime ([1]), the observed current 
is in fact an average (. . .)\^ over lengths d + where x 
ranges over Ac- This averaging effectively yields 



(1 + cos(2fci.d))A, = 1 



(17) 



so that the dependence on hp disappears. The results 
for the current presented henceforth are thus meant upon 
performing this averaging procedure. 
Secondly, the electron-electron interaction strongly af- 
fects the behavior of the correlators C^^a" Although 
the full expression for the latter is quite lengthy (see the 
appendix for details), important insights can be gained 
from the analysis of the scaling dimensions of the two 
operators; one obtains that for kpd ^ 1 



\C^!a'ikFd^O\ 



2-Sp 



(18) 



with(5p = (if+l/if + 2)/4and(5A = (l/if + 3)/4. While 
for vanishing interaction (K = 1) the scaling dimensions 
of the two processes coincide, the electron interaction 
modifies the power laws of these two processes in a dif- 
ferent way: the contribution of A processes decays faster 
than the one of P processes {5p < 6a)- Remarkably, 
this does not imply that for a sufficiently long junction 
the Josephson current is mainly due to P processes. 
Indeed, an electron pair traveling along the junction also 
acquires a phase, which results into oscillating factors 
in the correlators. Since the Josephson current (fT4|) 
depends on the integral over the imaginary time variable 

P j A 

^, not only the decay rate but also the phase of C^^^t' 
matters. Since the dynamics of the electrons is coupled 
by the interaction, these phase factor are also affected 
by the value of K. 

In the case that electron interaction is neglected (K — 
1), the effect of alternating phases is so strong that the 
total contribution of P processes vanishes. Indeed, when 
integrating over all possible pair momenta, the phase ac- 
quired by the electron pair traveling along the junction 
oscillates, yielding a cancellation of the different contri- 
butions, except for those processes in which the total 
pair momentum is vanishing. While this condition can 
be fulfilled by A processes, simple phase-space consider- 
ations show that the total contribution of P processes 
is suppressed. As a consequence, the Josephson current 
through a SWNT is predicted to be independent of the 
filling momentum qp^ and one obtains 




10- 



27re 



3c 



sinh(27re) 



(19) 



Figure 3: (solid curve): The oscillations of the dimension- 
less Josephson current as a function of junction length d for 
a S-I-SWNT-I-S junction, at zero temperature. The SWNT 
filling factor is qp/kp = 10~^7r/8 and its interaction parame- 
ter is A' = 0.3. The oscillations have a period 2qFd and decay 
with a power- law with an interaction-dependent exponent 5p. 
(dashed curve): The contribution of A processes to the crit- 
ical current is monotonous and positive, indicating that the 
oscillations originate from P processes (see text). 



At zero temperature, jc = 1, and the Josephson current 
scales as \/d due to the /o term (see Eq. ([13])), whereas 
at finite temperature it is exponentially suppressed. 

In contrast, when electron-electron interaction is 
taken into account {K ~ 0.3), the dynamics of the 
two electrons forming any pair is coupled, and the 
mechanism leading to the cancellation occurring in the 
non-interacting case is not valid. Electronic correlations 
both affect the contribution of A processes and make an 
oscillating contribution in 2qpd arise from P processes. 
These oscillations are characterized by a much longer 
period than the one discussed above, and may become 
observable if ^i^Ac ^ 1, a condition which is definitely 
realistic: The value of the filling factor qp jkp can indeed 
be adjusted by an external gate bias, and the recent 
developments in contact technology allow to realize 
extremely thin contacts, such as finger-shaped electrodes 
of about lOnm, or superconducting STM tips. Here we 
show that in this case the Josephson current exhibits 
interesting novel features. 

Fig. [3] displays the dimensionless Josephson critical 
current jc, Eq. (|14p . as a function of the junction length 
d for a SWNT with interaction strength K = 0.3 at 
zero temperature. We recall that the approximation of 
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non-interacting electrons would predict a constant value 
for jc- In contrast, in a SWNT the strong electron inter- 
action leads jc to decay with an oscillatory behavior as a 
function of d. While the power law decay has been pre- 
dicted also for usual two-channel LL, the oscillatory be- 
havior is purely due to the four channel band structure of 
nanotube. Importantly, this implies that the sign of the 
Josephson current depends on the length of the junction, 
and that SWNT can be used to reahze a 7r-state. We 
emphasize that this effect originates from P processes; 
the contribution to jc due to the A processes, described 
by the dashed curve in Fig. [31 is indeed monotonous and 
always positive. 

Fig. m shows jc for a S-I-SWNT-I-S junction with length 
d = 6 • lO^kp^ - 360 nm: the Josephson current oscil- 
lates with a period Tr/kpd as a function of the filling 
factor qp/kp, around a value (dashed line) which repre- 
sents the contribution of ^ processes, independent of qp. 
By tuning qp with a gate voltage, the switching from a 
to a TT-junction behavior can be induced. 
Finally, Fig. [5] shows the behavior of function 
of the dimensionless temperature O = ksT d/hvp. As 
expected from Eq. (|14p . thermal fluctuations suppress 
the Josephson effect at a temperature of the order of 
ksT ~ hvp/d (for a 100 nm long junction this corre- 
sponds to r = 60°K); the figure elucidates the crucial 
role played by the interaction in determining the relative 
magnitude of P processes with respect to A processes: 
while for a non-interacting wire {K = 1) the P processes 
contribution vanishes, for a SWNT {K ~ 0.3) the latter 
dominate for a sufficiently long junction. 



IV. CONCLUSIONS 

In this paper we have investigated the dc Josephson 
effect in a S-I-SWNT-I-S junction. The effects of the 
electron-electron interaction on the critical current jc 
have been particularly analyzed by using the Luttinger 
Liquid theory. We have found that jc oscillates with a 
factor 2qpd, where d is the junction length and qp the 
Fermi momentum characterizing the electron filling with 
respect to the band crossing point kp. These oscilla- 
tions are a signature of the peculiar band structure and 
of the strong electronic correlations present in SWNT. 
We emphasize that they would indeed not appear in 
non-interacting systems. Remarkably, this effect implies 
that ballistic SWNT can be used to realize tunable tt- 
junctions, for the sign of the critical current can be con- 
trolled by varying either the filling factor or the junction 
length (see Figs. [3] and |4]). The former can be tuned 
through an external gate voltage. The latter can be 
changed for instance by moving the superconducting tip 
of an STM2^i^i2i along the nanotube. The typical value 
of kp is of the order of 20 nm~^, so that the predicted 
oscillations should be observable in junctions with length 
d > lOOnm operating at temperature of the order of some 
°K or below; the SWNT should have an electron filling 
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.4 I I I < I < I I 

2 4 6 

Figure 4; (solid): The dimensionless critical current jc as a 
function of the filling factor qp/kp at zero temperature for 
a S-I-SWNT-I-S junction with length d = 6 • lO^fc^V The 
interaction parameter is K = 0.3. By tuning the electron fill- 
ing, e.g. with a gate voltage, the sign of the critical current 
can be reversed, tuning the junction from a into a vr-state. 
(dashed) : The contribution of ^-processes to the critical cur- 
rent, independent of the filling factor. 



momentum qp ranging from up to a small fraction of 

Xc\ 



Acknowledgements 

Fruitful discussions with F. Giazotto as well as finan- 
cial support from HYSWITCH EU Project are greatly 
acknowledged. 



Appendix A: Ol CORRELATION FUNCTIONS 

In this appendix we provide the expressions for the 
T-ordered correlation functions appearing in the com- 
putation of the Josephson current Eq. HH). The corre- 
lation functions can be written as the product of their 
ground state value and a thermal fluctuations contribu- 
tion, which equals 1 at zero temperature. Explicitly: 

c^{kpd, = [q^ («^) {kpd, (^^) {kpd, 0] ' 

(Al) 

{kpd, = Isf («^) {kpd, flf ^^^^ {kpd, 1 ' 
where: 
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a 








a + z'l^ 






-zrJ 



(A3) 



0.5 



© 



1.5 



with z"^ = rsign(^) + l^l/if) and d = kpa. The 

contribution due to thermal fluctuations reads 



K 



sin(if7rez^) 

ttQzI* 
sin(7re^r) 



Sa-I/2 



K 



sin(/s:7rez^) 



(A4) 



Figure 5: (solid): The dimensionless critical current jc as a 
function of the reduced temperature for a junction with 

length d = 6 • lO^fc^S rcahzcd with a SWNT with filling fac- 
tor qp/kp ~ 10^'^7r/12 and interaction parameter K = 0.3. 
(dashed): The contribution of A-processcs to the critical cur- 
rent. Differently from the case of a non-interacting wire, in a 
SWNT the contribution due to tunnel of pairs with antipaxal- 
lel momenta {A processes) is dominated by the one originating 
from pairs with parallel momenta (P processes). 



A{GS) 



{kFd,0 



Sa-1/2 



a + 

d + z^ 
a + zH 



d + z^ 



0f(^^)(fc^d,O 



KirOzl; 



sin(i^7rez^) 

ttOzI* 
sin(7rez[*) 



K-rrQz'ii 
sin(^7rG^ 



The cut-ofl' d rcnormalizcs the bare tunneling amplitude 
T in different ways for P and A processes. In particu- 
lar, one has r ^ % = T{kFa)'^^^*-'^\ with i = A,P. 
(A2) Typically a < 1 (here we have chosen a = 0.5). 
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